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<^ ■ Abstract. An attempt is made to guess the overall cosmic abundance 

^^ \ of 'metals' and the contribution made by the energy released in their 

production to the total intensity of extragalactic background light (EBL). 
With a comparable or somewhat larger amount coming from white dwarfs, 
and a probably quite modest contribution from AGNs, one can fairly 
easily account for the lower end of the range of existing estimates for 
the total EBL intensity (50 to 60 nwt m -2 sterad -1 ), but it seems more 
difficult should some higher estimates (90 to 100 in the same units) prove 
to be correct. 
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t^" ' 1. Introduction 

O" 

There are certain more or less well or badly determined integral constraints on 
the past history of star formation in the universe. These include 

»^-( • Cosmic baryon density S7b, now fairly well determined both from pri- 

mordial deuterium (O'Meara et al 2000) and from the CMB fluctuation 
£j ', spectrum (e.g. Turner 2001). 

• Cosmic mass density of stars Q*, rather less well determined as it involves 
a combination of luminosity-density measurements with an IMF and evo- 

k> i lutionary population synthesis models. 

C^ ■ • Extragalactic background light (EBL) intensity, now known within a factor 

of 2 or so from COBE (FIRAS and DIRBE) and galaxy counts in the 
optical and near IR (e.g. Gispert, Lagache & Puget 2000). 

• Cosmic abundance of 'metals', fiz, due to the heavy-element content of 
stars, the interstellar medium and the intergalactic medium, a quantity 
that is very poorly known and largely a matter of guesswork. In this talk 
I shall nevertheless make some guesses, so that at least one can see more 
easily how things relate to one another. In particular, the metallicity 
Zigm of the intergalactic medium has tended to be either neglected or 
underestimated in models of the past star formation rate, and it is of 
interest to ask about its relation with EBL. 
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2. The cosmic inventory 

A useful starting point is the cosmic baryon budget drawn up by Fukugita, 
Hogan & Peebles (1998), hereinafter FHP, shown in the accompanying table. 
The total from Big Bang nucleosynthesis (BBNS) adopted here agrees quite 
well with the amount of intergalactic gas at a red-shift of 2 to 3 deduced from 
the Lyman forest, but exceeds the present-day stellar (plus cold gas) density by 
an order of magnitude. [] 

FHP pointed out that a dominant and uncertain contribution to the baryon 
budget comes from intergalactic ionized gas, not readily detectable because of 
its high temperature and low density. The number which I quote is based on 
the assumption that the spheroid star-to-gravitational mass ratio and baryon 
fraction are the same in clusters and the field, an assumption that had also 
been used previously by Mushotzky & Loewenstein (1997). The resulting total 
star-plus- gas density is within spitting distance of Ob from BBNS, but leaves 
a significant-looking shortfall which may be made up by some combination of 
MACHOs and low surface-brightness galaxies; it is not clear that a significant 
contribution from the latter has been ruled out (cf O'Neil 2000). 

3. Global abundances and yields 

We now have the tricky task of estimating the total heavy-element content of 
the universe. Considering stars alone, it seems reasonable to adopt solar Z as 
an average, but the total may be dominated by the still unseen intergalactic 
gas, which Mushotzky & Loewenstein argue to have the same composition as 
the hotter, denser gas seen in clusters of galaxies, i.e. about 1/3 solar .^] It could 
be the case, though, that the metallicity of the IGM is substantially lower in 
light of the metallicity-density relation predicted by Cen & Ostriker (1999) and 
in that of the low metallicities found in low red-shift Ly-a clouds by Shull et al 
(1998). Against this, we have neglected any metals contained in LSB galaxies 
or whatever makes up the shortfall between Otgm and Q,b, so we are being 
conservative in our estimate of Viz- 

The mass of heavy elements in the universe is related to that of stars through 
the yield, defined as the mass of 'metals' synthesised and ejected by a generation 
of stars divided by the mass left in form of long-lived stars or compact remnants 
(Searle &; Sargent 1972). The yield may be predicted by a combination of an 
IMF with models of stellar yields as a function of mass, or deduced empirically 
by applying a galactic chemical evolution (GCE) model to a particular region 



lr The stellar density taken here from FHP is based on B-luminosity density estimates and might 
be revised upwards by 50 per cent in light of SDSS commissioning data (Blanton et al 2000) 
or downwards by 20 per cent in light of 2dF red-shifts plus 2MASS K-magnitudes (Cole et al 
2000); in either case we are following FHP in assuming the IMF by Gould, Bahcall & Flynn 
(1996), which has 0.7 times the M/Lv ratio for old stellar populations compared to a Salpeter 
function with lower cutoff at O.15M0. 

2 This refers to iron abundance, the relation of which to the more energetically relevant quantity 
Z is open to some doubt. Papers given at this conference indicate an SNIa-type mixture 
in the immediate surroundings of cD galaxies with maybe a more SNII-like mixture in the 
intra-cluster medium in general; for simplicity I assume the mixture to be solar. 
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Table 1. Inventory of cosmic baryons and 'metals' 

Densities expressed as Q, in units of p cvl t = 1.54 x 10 n /iyo Mq Mpc 

All baryons from BBNS 

(D/H = 3 x 1(T 5 a ) 0.04 /i^ 2 

Stars in spheroids 0.0026 h^ 

Stars in disks 0.0009 h^ b 

Total stars 0.0035 h^ b 

Cluster hot gas 0.0026 /if L5 b 

Group/field hot gas 0.014 /if L5 b (0.004/if 5 1 in O vi systems c ) 

Total stars + gas 0.021 /if 1 ' 5 b 



Machos + LSB gals 



'70 



Q z (stars, Z = 0.02 d ) 7 x 10~ 5 h^ 

n z (hot gas, Z = .006) 1.0 x 10~ 4 /if L5 b 

1.2 x 10~ 4 /ifo 1 - 3 e 

Yield pz/p* 0.051 /ifp - 3 (~ 3Z Q \ 



Damped Ly-a (H i) 0.0015 h^ b 'f 

Ly-q forest (H+) 0.04 h^ b '9 



Gals + DM halos 

(M/L = 210/i 70 ) 0.25 6,/l 



All matter 



i.5\ n q^-O-S 6,-, 



(/B = .056/ 1 " i& ) 0.37/i. 



70 



a O'Meara et al 2001; but see also Pettini & Bowen 2001; b Fukugita, Hogan & Peebles 
1998; c Tripp, Savage & Jenkins 2000; d Edmunds & Phillipps 1997; e Mushotzky & 
Loewenstein 1997; ? Storrie-Lombardi, Irwin & MacMahon 1996; 9 Rauch, Miralda- 
Escude, Sargent et al 1998; h Bahcall, Lubin & Dorman 1995; l White & Fabian 1995. 
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Figure 1. Oxygen abundance distribution function in the solar neigh- 
bourhood, after Pagel & Tautvaisiene (1995). 

like the solar neighbourhood and comparing with abundance data. E.g. Fig 1 
shows an abundance distribution function for the solar neighbourhood plotted in 
a form where in generic GCE models the maximum of the curve gives the yield 
directly, and it is a bit below Z Q . Similar values are predicted theoretically using 
fairly steep IMFs like that of Scalo (1986). In Table 1, on the other hand, if 
we divide the mass of metals by the mass of stars, we get a substantially higher 
value, corresponding to a more top-heavy IMF. 

There are two other indications for a top-heavy IMF, one local and one 
in clusters of galaxies themselves. The local one is an investigation by Scalo 
(1998) of open clusters in the Milky Way and the LMC, where he plots the IMF 
slopes found as a function of stellar mass. The scatter is large, but on average 
he finds a Salpeter slope above O.7M and a virtually flat relation (in the sense 
dN/dlogm ~ 0) below, which could quite easily account for the sort of yield 
found in Table 1. The other indication is just the converse of the argument we 
have already used in guessing the abundance in the IGM: the mass of iron in 
the intra-cluster gas is found (Arnaud et al 1992) to be 



M Fe /L v = O.O18M /L G 



0> 



1) 



where Ly is the luminosity of E and SO galaxies in the cluster. As has been 
pointed out by Renzini et al (1993) and Pagel (1997), given a massilight ratio 
less than 10, we then have 



M Fc (gas)/M* > 1.8 x 10~ 3 = 1.6Z (Fe) 
M Fe (*)/M* ~ Z (Fe) 
Yield = M Fe (total) /M* > 2.6Z Q . 



(2) 
(3) 
(4) 



The argument is very simple; the issue is just whether such high yields are 
universal or confined to elliptical galaxies in clusters. 

4. Cosmic star formation and chemical evolution: GCE vs HDF 



The deduction of past star formation rates from rest-frame UV radiation in 
the Hubble and other deep fields as a function of red-shift is tied to 'metal' 
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production through the Lilly-Cowie theorem (Lilly & Cowie 1987): 

1 



PlO) 



vkPv\jv(z) = -L Fm (z) 



0.007c 2 pz(z)(l + a) (5 
0.018c 2 p z (conv.), 



-i 



(5) 

(6) 
(7) 



where (1 + a) ~ 2.6 is a correction factor to allow for production of helium as 
well as conventional metals and f3 (probably between about 1/2 and 1) allows for 
nucleosynthesis products falling back into black-hole remnants from the higher- 
mass stars, e is the fraction of total energy output absorbed and re-radiated by 
dust and fn is the frequency at the Lyman limit (assuming a flat spectrum at 
lower frequencies). The advantage of this formulation is that the relationship is 
fairly insensitive to details of the IMF. 
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Figure 2. Global comoving star formation rate density vs. lookback 
time compiled from wide-angle ground-based surveys (Steidel et al. 
1999 and references therein) assuming E-de S cosmology with h = 0.5, 
after Pettini (1999). Courtesy Max Pettini. 

Eq (7) is the same as eq (13) of Madau et al (1996), so I refer to the 
metal-growth rate derived in this way as p\ (conventional). 

Assuming a Salpeter IMF from 0.1 to IOOMq with all stars above IOMq ex- 
pelling their synthetic products in SN explosions, one then derives a conventional 
SFR density through multiplication with the magic number 42: 



p*(conv.) = 42p z (conv.) = 42p L /.018c 2 . 



In general, we shall have 



p*(true) = 7 / o*(conv.), 



(8) 



(9) 



where 7 is some factor. E.g., for the IMF adopted by FHP, 7 = 0.67, whereas 
for the Kroupa-Scalo one (Kroupa et al 1993) 7 = 2.5. 
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Finally, the present stellar density is derived by integrating over the past 
SFR and allowing for stellar mass loss in the meantime, and the metal density 
is related to this through the yield, p: 



p*(0, true) 
/9 Z (0,true) 



«7 / p*(conv.)dt 



p/9*(0,true) 



2.6 
1 + a 



(10) 

(3 r 
— J p*(conv.)dt (11) 



(where a is the lockup fraction), whence (if a = 1.6) 



P 



P 



42«7' 



(12) 



which can be compared with Zq ~ 1/60. It was pointed out by Madau et al 
(1996) that the Salpeter slope gives a better fit to the present-day stellar density 
than one gets from the steeper one — a result that is virtually independent of 
the low-mass cutoff if one assumes a power-law IMF. 

Eq (8), duly corrected for absorption, forms the basis for numerous discus- 
sions of the cosmic past star- formation rate or 'Madau plot'. Among the more 
plausible ones are those given by Pettini (1999) shown in Fig 2 and by Rowan- 
Robinson (2000), which leads to similar results and is shown to explain the far 
IR data. Taking 7 = 0.62 (corresponding to a Salpeter IMF that is flat below 
0.7 Mq) rather than Pettini's value of 0.4 (for an IMF truncated at 1M ), and 
a = 0.7, we get the data in the following table. 



Table 2. Inventory of stars and metals at z = and z = 2.5 





z = 


z = 2.5 


P* = ot 7 / p* (conv. ) dt 


3.6 x 10 8 M Mpc~ 3 


9 x 10 7 M Mpc~ 3 


n* = p*/1.54 x io n /4 
0,(FHP 98) 


.0024/if 2 
.0035/ifo 1 


6 x lO^/ifo 2 


Pz = PP* = /3p*/(42«7) 
Vl z (predicted) 


2.0 x lO 7 /3M Mp C - 3 
1.3 x lO- 4 /^ 2 


5 x 10 6 /3M©Mpc" 3 
3.2 x lO- 5 /^ 2 


Viz (stars, Z = Z Q ) 
Viz (hot gas, Z = 0.3Z Q ) 


7 x lcr 5 /^ 1 

1.0 x 10- 4 /if L5 
^0.5</3< 1.3 




Vl z (DLA, Z = 0.07Z Q ) 

Vl z (Ly. forest, Z = O.OO3Z ) 

Viz (Ly- break gals, Z = 0.3Z Q ) 

Viz (hot gas) 




2 x lO^/ifo 1 

1 x 10- 6 /if 2 
? 
? 
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Table 2 indicates that the known stars are roughly accounted for by the 
history shown in Fig 2 (or by Rowan- Robinson) and the metals also if (5 is close 
to unity, i.e. the full range of stellar masses expel their nucleosynthesis products. 
At the very least, (3 has to be 1/2, to account for metals in stars alone. The 
other point arising from the table, made by Pettini, is that at a red-shift of 2.5, 
1/4 of the stars and metals have already been formed, but we do not know where 
the resulting metals reside. 

5. Extragalactic background light 
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Figure 3. Spectrum of extragalactic background light, based on 
COBE data after Hauser 2001 (diamonds with error bars, dotted and 
short-dash curves), Madau & Pozzetti 2000 (squares), Totani et al 2001 
(crosses), Bernstein, Freedman & Madore 2001 (triangles) and Armand 
et al 1994 (asterisk). The broken-line curve (Biller et al 1998) and hor- 
izontal dash-dot line (Hauser 2001) show upper limits based on lack 
of attenuation of high-energy 7-rays from AGNs and the solid curve 
is from the model by Pei, Fall & Hauser (1999). The arrow showing 
an upper limit at 10[im is from unpublished thesis work by A. Barrau, 
cited by Gispard, Lagache & Puget (2000). 

Fig 3 shows the spectrum of extragalactic background light (EBL) with the 
model fit by Pei, Fall & Hauser (1999). Gispert, Lagache & Puget (2000) have 
estimated the total EBL / I v dv based on observation to lie within the following 
limits: 

A < 6/um: 20 to 40 nwt m~ 2 sterad -1 
A > 6/zm: 40 to 50 " " 
Total: 60 to 90 " " 
(The total from the model of Pei, Fall & Hauser (1999) is 55 in these units.) 
We use the estimates of stellar and metal densities in Tables 1 and 2 together 
with eq (7) and an assumption about the mean red-shift of metal formation to 
derive the EBL contributions from: 



• Metals in stars; Z = Z & = 0.02: 

Pz(*) = 7 x 10~ 34 h 70 gm cm" 



(13) 
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, = »£f = ( ^ rl 3xl034 

47r(l + z) 1 + z 

~ 9/3 -1 /170 nwt m~ 2 sterad _1 . (15) 

• Metals in diffuse gas/dark baryons, ^ ~ 0.007: 

/5z = ^zPcrit = 1.2xKT 33 /i 7 o - 7 /? gmcm- 3 . (16) 

I = 16/i?5 7 nwtm _2 sterad _1 . (17) 

• Helium, carbon etc in white dwarfs and red-giant interiors. 

Here we use eq (6) without the (1 + a), since most of the nuclear energy 
is already released on reaching this stage. Assuming most stars to belong 
to an old population so that 

PWD ^ O.lp* = 3.5 x 10 33 /i 70 gmcm -3 , (18) 

0.007c 3 p W D , 3 33 

1 = a /i i \ = T^~ 5x 10 IWD- 19 

47r(l + z) 1 + z 

= 18/i 70 nwt m~ 2 sterad _1 . (20) 

• AGN contribution 

Madau & Pozzetti (2000) and Brusa, Comastri & Vignali (2001) have made 
estimates of the AGN contribution to EBL based on the the abundance 
of massive black holes and that of obscured hard X-ray sources, respec- 
tively. They agree that the contribution is quite small, of order 5 nwt m~ 2 
sterad -1 . 

The upshot is that these readily identifiable contributions add up to 48 nwt 
m -2 sterad -1 , well within range (given the obvious uncertainties in mean z and 
other parameters) of the lower estimate given at the beginning of this section. 
It is interesting to note that white dwarfs and intergalactic metals come out 
as the major contributors, either one predominating over metallicity in known 
stars. To reach the higher estimate may involve some more stretching of the 
parameters. 

I thank Richard Bower, Jon Loveday and Max Pettini for helpful informa- 
tion and comments. 
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